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Abstract 

The definition of order indices for density matrices is extended to finite systems. This 
makes it possible to characterize the level of ordering in such finite systems as macro- 
molecules, nanoclusters, quantum dots, or trapped atoms. The general theory is exempli- 
fied by explicit calculations of the order index for the first-order density matrix of bosonic 
atoms confined in a finite box at zero temperature. 
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1 Introduction 



We collaborated with John A. Coleman for 20 years, since 1990 till his death in 2010. Moreover, 
he was not merely our co-author, but during this time, he became our very close friend, with 
whom we spent many hours discussing various problems. Wc are keeping warmest memories of 
John and his wife Marie Jeanne. With great respect, we devote this paper to the memory of our 
friend and colleague John Coleman. 

One of the interesting ideas, we developed with John Coleman, was the notion of the order 
indices for density matrices, introduced in Ref . [1] . The properties and applications of the order 
indices to different types of bulk matter, considered in thermodynamic limit, were studied in 
Refs. [2-5] and summarized in the book [6] . For infinite systems, however, it is possible to define 
the standard order parameters (see, e.g., [7,8]), because of which the use of the additional notion 
of the order indices could seem to be unnecessary. Nevertheless, as has been shown [1-3], the 
order indices are useful even for infinite systems, where a kind of mid-range order arises. 

The principal difference of the order indices from the order parameters is that the former 
can be introduced not only for infinite systems, but also for finite systems. In recent years, the 
investigation of finite systems has become of high importance due to the widespread technological 
applications of various finite objects. As examples, we can mention quantum dots, metallic grains, 
different granular materials, nanoclusters, trapped atoms, and a variety of macromolecules, in- 
chiding biomolecules. For finite systems, as is well known, the order parameters are not defined 
[7,8]. In that case, the order indices can become principally important, as far as they can be 
defined for systems of any size. Then the order indices could characterize the amount of or- 
der specific for finite systems. It is the aim of the present paper to extend the definition and 
application of the order indices for finite systems. 

In Sec. 2, we introduce the notion of order indices for arbitrary operators, which is specified 
for generalized density matrices in Sec. 3. The application to the usual reduced density matrices 
of statistical systems is given in Sec. 4, without invoking thermodynamic limit. In Sec. 5, we 
exemplify the consideration for bosonic atoms trapped in a finite box. Calculations for the order 
index of the first-order density matrix are given in Sec. 6, where the order index is treated as a 
function of the number of particles and of atomic interactions. Section 7 concludes. 



2 Operator order indices 

Order indices can be introduced for operators of any nature [9] . Let A be an operator acting on 
a Hilbert space H. The operator is assumed to possess a norm and a trace TrA, with the 
trace taken over the space H. In all other aspects, it can be arbitrary. The operator order index 
is 

log|Tr^| 

The logarithm can be taken with respect to any base, since log^ z — log^ log^ a. This definition 
connects the norm and trace of the operator through the relation 

jjijl = |Tri|'^(^^ (2) 

An operator Ai is said to be better ordered than A2, if and only if 

> uiA^) . (3) 
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Respectively, two operators, Ai and A2 are equally ordered, provided that uj{Ai) = u{A2). 

The operator norm can be defined as a norm associated with the vector norm \ip\ for a nonzero 
vector if eT-L, so that 

\A^\ 



\A\ I = sup 



(4) 



Employing the scalar product {ip, if) for defining the norm yields the Hermitian vector norm 



(fi, (fi). The corresponding Hermitian operator norm is 



\A\ 



sup 



'iAip,Aipy 


1/2 




1/2 




= sup 







For an orthonormalized basis {(pk}, labelled with an index k, in T-L, such that 

n = Spanjt{|(^fe)} , 

the Hermitian norm (5) becomes 

1/2 



sup 

k 



{A(pk:A(pk) 



(5) 



(6) 



(7) 



If the operator A is self-adjoint, then its Hermitian norm simplifies to 



sup 



sup \{(pk:A(pk)\ 
k 



(8) 



The eigenf unctions of a self- adjoint operator, defined by the eigenproblem 

Aipk = Akcpk , 

form an orthogonal basis that can be normalized. The space basis can be chosen as the set of 
these eigenfunctions of the considered operator. Then the Hermitian norm becomes the spectral 
norm 



\\A\ \ — sup \ Ak\ 

k 

When the operator is semi-positive, then \Ak\ = A^. 
For a semi-positive operator. 



(9) 



Therefore, for such an operator. 



|i||<Tri (i>0). 



u{A) < 1 (i > 0) . 



(10) 



In that way, the order index (1) makes it possible to characterize the level of order in operators 
and to compare the operators as being more or less ordered. Instead of the Hermitian norm, one 
could employ some other types of operator norms, for instance, the trace norm [9]. But the use 
of the Hermitian norm is more convenient for physical and chemical applications. 
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3 Generalized density matrices 

Self-adjoint operators play a special role in applications, defining the operators of observable 
quantities. One can resort to the coordinate representation, defining the physical coordinates 
through X, implying the set of all variables characterizing a particle. These can include Cartesian 
coordinates, spin, isospin, component-enumerating labels, and like that. The arithmetic space of 
all admissible values of the physical coordinates is denoted as A* = {x}. 

Let A{x) be an operator of a local observable from the algebra of local observables A = {A{x)} 
given on the Fock space J^. And let = be the algebra of field operators on 

describing the system. The direct sum of these algebras is an extended local algebra 

Aext = A A^ . ( 1 1 ) 

The set x"' = {xi, X2, ■ ■ ■ , x„} of the coordinates of n particles pertains to the arithmetic space 

that is an n-fold tensor product. The differential measure on the above space X"- is defined as 

dx"' = dxidx2 ■ ■ ■ dxn ■ 
A function ip{x^) can be treated as a vector 

= [(^(x")] e Un (12) 
in a Hilbert space l-Lm where the scalar product is given by 

^l^r.^ U*{x'^)^{x-)dx- . (13) 



We consider a quantum system, whose state is given by a statistical operator p being a 
semi-positive self-adjoint operator normalized as 

TV^p =1 (p+ = p > 0) . (14) 

Taking any representative A[x) of the extended algebra (11), we can define a matrix 

Dl=[DA{x^,y^)] , (15) 

which is a matrix with respect to the variables x, with the components 

DA{x\y^)=TvrA{x^)...A{x^)pA+{yr,)...A+{y^) . (16) 

The action of matrix (15) on vector (12) is defined as the vector with the components 



j DA{x\y^^)^{y^)dy^^ 



(17) 



As is seen from construction, matrix (15) is self-adjoint and semi-positive, because of which it 
can be termed the generalized density matrix. 
The norm of matrix (15) is defined as 

irai=sup^i|^. (18) 



And the trace of this matrix is 

TrlDl = J dx" . (19) 

The order index of the generahzed density matrix (15) is given by the expression 

uiPn - ^"gll^^ll (20) 
log|TrL)^| ' ^ ^ 

with the norm and trace defined as above. 

4 Reduced density matrices 

A particular case of the generalized density matrices, being the most important for applications, 
is that of reduced density matrices. An n-th order reduced density matrix is the matrix 

p„ = [p(x",y")], (21) 

with the components 

p(x",y") = Tr^V'(a:i) . . . V^(a:„)pV't(^^) . . .^t(^^) . (22) 

The relation of the reduced density matrix with the generalized density matrices is given by the 
equations 

Similarly to definition (20), the order index of the reduced density matrix (21) is 

N log||pn|| 

In this definition, we do not take thermodynamic limit, as in Refs. [1-3]. 
The eigenproblem 

Pn^nk = Nnkfnk (24) 

yields the eigenvalues 

Nnk = VnkPnVnk (25) 

that define the spectral norm 

llPnII = supA^„fc . (26) 

k 

And the trace of the n-th order density matrix is given [6] by the normalization 

We keep in mind a finite system with a finite number of particles N. This number is assumed 
to be sufficiently large, N ^ 1, but finite. For a large N and fixed n <^ N, we can invoke the 
Stirling formula, yielding 

Trp„ ~ (^-j . (28) 
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Taking the natural logarithm gives 

lnTrp„ = n(lnA^-l) . 

Then the order index (23) reads as 

From the properties of the reduced density matrices [6], we know that, in the case of Bose 
particles, 

\\Pn\\<{bnNr, (30) 

where 6„ is a constant. This results in the inequality 

^(M < '-^^^^ (Bose) . (31) 
For Fermi particles, depending on the order of the density matrices, we have 

||p2n-l|| < (c2n-liV)-l , ||p2n|| < (c2nA^)" , (32) 



because of which 



, (n-l)(ln7V + lnc2n-i) 

^^^^'^-^^ - (2n-l)(lnAr-l) 



lnAr + lnc2n .„„^ 
^^^^"^ 2(lnAr- 1) (^ermz) . (33) 

Taking into account large N results in the inequalities for Bose particles, 

uo{pn) < 1 (Bose) , (34) 

and for Fermi particles, 

n — 1 1 

^(P2n-i) < 7, 7 ' ^{p2n) < 7: {Fermi) . (35) 

Zn — L 2 

5 Bose system with quasi-condensate 

To give a feeling how the order indices describe physical ordering in finite systems, let us consider 
a cloud of N Bose atoms trapped in a finite box of volume V. Systems of trapped Bose atoms are 
nowadays intensively studied both theoretically as well as experimentally [10-24]. We consider 

the case of either spinless atoms or that where atomic spins are frozen by an external magnetic 
field, so that the spin degrees of freedom are not important, but only the spatial variables r are 
considered. 

As is known, in an infinite Bose system at low temperature, under thermodynamic limit, when 
— )■ oo, there appears Bose- Einstein condensate. But if the system is finite, no matter how large 

it is, there can be no well defined phase transition, hence, no Bose-Einstein condensation [20,24]. 
In a finite system, there can exist only a kind of quasi-condensatc [25]. Below, we demonstrate 
the calculation of the order index for a Bose system in a finite box, at low temperature, when a 
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quasi-condensate appears. For brevity, we shall often use the term condensate, keeping in mind 
quasi-condensate . 

We start with the standard energy Hamiltonian 



J ^\r) ^(r) dr + ij\r)i!\r')<^{r - r')^(r')^(r) drdr' , 



(36) 



expressed through the field operators ■0(r) satisfying the Bose commutation relations. Here and 
in what follows, the system of units is used, where the Planck and Boltzmann constants are set 

to one, h = 1, ks = 1- 

To take into account the appearance of quasi-condensate, we employ the Bogolubov shift [26] 
of the field operator 

V-W^TyW + ViW, (37) 

separating it into the condensate function rj and the operator of the normal, uncondensed, par- 
ticles ipi. To avoid double counting, the condensate and normal degrees of freedom are assumed 
to be orthogonal, 

y^77*(r)V'i(r) dr = . (38) 
The operator of uncondensed particles satisfies the conservation law 

(^i(r)) = 0. (39) 
There are two normalization conditions, for the number of condensed particles, 

No = J \v{r)\'dr, (40) 

and for the number of uncondensed particles, 

N, - (A) , (41) 
with the number operator of uncondensed atoms 

N, = J 4(r)M^)dr . (42) 

The total number of atoms in the box is the sum 

N^No + Ni. (43) 

To guarantee the validity of the normalization conditions (40) and (41), as well as the con- 
servation law (39), it is necessary to introduce the grand Hamiltonian 

H = H- iioNo - iiiNi - A , (44) 

in which 

A = y [X{r)4{r) + A*(r)^i(r)] dr , (45) 

with /Iq, fJ,i, and A(r) being the Lagrange multipliers. 

The equations of motion are given by the variational equations for the condensate function 



and for the field operator of uncondensed atoms 

where t is time and the angle brackets imply statistical averaging. These equations are equivalent 
to the Heisenberg equations of motion [24,27]. 
The first-order reduced density matrix is 

p(r,r') = (V^^(r')V^(r)). (48) 

This, in view of the Bogolubov shift (37), reads as 

p(r, r') = v*ir')vir) + (V^I (r')^i(r)) . (49) 

The cigenfunctions of the reduced density matrix (48) are called natural orbitals [6] . If these 
eigcnfunctions are denoted as v^jfc(r), with a labelling quantum multi- index k, then the spectrum 
of the density matrix is defined by the quantities 

N,k = / ^l{r)p{r, r')<^fc(r') drdr' . (50) 

According to expression (49), there are two terms in the latter integral. One term. 



(51) 



V*{r)'Pk{r) dr 

characterizes the condensed atoms, while another term, 

n, = J ^l(v){^|;l(r')Mr))M^') dvdr' , (52) 
defines the distribution of uncondensed atoms. With the notation 

«fe = / <^fe(r)V'i(r) dr , (53) 

this distribution takes the form 

Uk = (4«fc) • (54) 
In this way, the spectral norm of the density matrix (48) can be represented as 

I I = sup Nik = sup(7Vfc + Uk) . (55) 

k k 

In a similar way, it is possible to find the norms of the higher-order density matrices. Thus, 
for the second-order density matrix, we would have to find out the pairon spectrum [28]. But 
here we concentrate on the properties of the first-order density matrix. 
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6 Order index behavior 



Let us study the behavior of the order index 

of the first-order density matrix (48). For atoms in a box of volume V, the natural orbitals are 
the plane waves 

labelled by the wave vector quantum number k. The condensate wave function reduces to a 
constant 77 = ^JNq/V . Then the condensate spectrum (51) is 

Nk = iVo4o • 

The matrix eigenvalues (50) read as 

Nik = SkoNo + (1 - 6ko)nk . (57) 

And norm (55) becomes 

IIpiII = sup{A?o, sup rife} . (58) 

k 

To accomplish explicit calculations, we need to fix the form of the interaction potential enter- 
ing Hamiltonian (36). We shall keep in mind dilute Bose gas for which the interaction potential 
is well modelled by the local form 

$(r) = $o5(r) , $0 = 47r ^ , (59) 

m 

where is scattering length. We shall accomplish calculations invoking the Hartree-Fock- 
Bogolubov approximation in the self-consistent approach using representative ensembles [23,24,29]. 
We introduce the notations 

iVo A^i 1 

k 

P = ^ = Po + Pi, <^i = ^^nk, (60) 

k 

defining the mean densities of condensed (po) and uncondensed (pi) atoms, and also the so-called 
anomalous average ai, whose modulus \ai\ describes the number of correlated atomic pairs. The 
corresponding atomic distributions, at temperature T, read as 

where 



ujk = mc^ + 7^ , £fe 
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The sound velocity s satisfies the equation 

mc^ = $o(po + ^i) • (62) 

Quasi-condensate can arise only at low temperature. For concreteness, we take zero temper- 
ature T = 0. Then, from the above formulas, we have 



TT 



(63) 



In the calculation of ai, wc employ dimensional rcgTilarization [24,29]. 

It is convenient to introduce dimensionless quantities simplifying the formulas. Atomic inter- 
actions are characterized by the gas parameter 

7 = ■ (64) 

Dimensionless sound velocity is 



mc 



(65) 



The condensate and anomalous fractions are 



no = - = -, a=-. (66) 

In this dimensionless notation, the equation for the sound velocity (62) reduces to 

= 47r7(no + a) , (67) 

the condensate fraction becomes 

no = l-^, (68) 
and the anomalous fraction is 

TT V 7r 

Numerical solution to these equations is shown in Fig. 1, where s,a, and no are presented as 
functions of the gas parameter (64). 

The distribution of the normal atoms Uk increases as k diminishes. The minimal value of k 
is prescribed by the box volume V, so that 

1 



This gives 

sup rife = J N^/^ . 

k ^ 



Therefore, norm (58) is represented as 

||pi||=sup{noAr, |iV^/'} . (71) 
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The order index (56), for an infinite system in the presence of condensate, in thermodynamic 
hmit — )■ oo, is exactly one, which corresponds to long-range order. But for a finite system, the 
order-index behavior is not so trivial, varying between zero and one, depending on the system 
parameters. 

In the case of asymptotically weak atomic interactions, when 7 — > 0, from Eqs. (67) to (69), 
it follows 

8 3/2 64 3 

8 0/0 32 o 

—= 7"^/^ H 7^^ , 

Vtt tt 

.c.20F7^/^ + f 7^ + ^7^/^ (72) 



Then the matrix norm is 
and the order index tends to 



IpiII ~ uqN (7^0), 
u[pi) ~ 1 



In TV 

Using expansions (72), we find 

As is seen, the order index reduces to one when either the interaction strength tends to zero or 
the system becomes infinite. But for a finite system of interacting atoms, the order index is less 
than one. 

In the opposite case of arbitrarily strong interactions, when 7 ^ cxd, we have 

^ -3 
no-^7 , 

(9^)V3 ^ ^ 3 1 / ^4X1/3 

4 ' 64 ' 128 V 3 y ' ' 

For a finite system, with a fixed number of atoms N, when the interaction strength exceeds the 
value 0.317Ar2/9, we get 

llplll-^ 

Therefore, the order index behaves as 



/-^ N 1 0-436 , , 
^(Pi) - 3 + (7 ^ 00) , (75) 

which is again less than one. 

Thus, we see that, for a finite system, the order index varies between unity, when the gas 
parameter tends to zero, and the limiting value (75), when the gas parameter tends to infinity. 
That is, a finite system can exhibit only a kind of quasi-long-range order, or mid-range order. 
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The overall behavior of the order index (56) as a function of the interaction strength 7 and the 
number of atoms in the box A^, is shown in Figs. 2 and 3. When tends to infinity, the order 
index tends to unity. In the scale of Fig. 3, this is not seen for 7 = 1, since uj{pi) for this 7 
decreases after N — 11. In order to demonstrate that this is just a temporary decrease, we show 
the behavior of the order index in a larger scale of N, in Fig. 4. As is seen, the order 

index diminishes till N = 200, after which it increases. This increase is rather slow. Thus, for 
= 10^, the order index reaches the value 0.9, and for N — 10^°, it becomes 0.98. So that it 
tends to unity as ^ 00. 

7 Conclusion 

The notion of order indices for reduced density matrices is extended to the case of finite systems. 
In such systems there can be no long-range order and the standard order parameters are not 
defined. Contrary to this, the order indices have the meaning for finite systems, characterizing 
the level of ordering in these objects. Examples of finite systems are quantum dots, metallic 
grains, nanoclusters, trapped atoms, and various macromolecules, including biomolecules. 

Generally, in finite systems, there can exist only a kind of mid-range order. The level of such 
an order is well described by order indices. The suggested approach is illustrated by calculating 
the order index for the first-order reduced density matrix of bosonic atoms at zero temperature, 
trapped in a finite box. The order index, depending on the values of the number of atoms and 
the strength of their interactions, varies between zero and one. 

The order indices can be defined for finite systems composed of several types of particles. 
For such a system of several components, it is possible to define the order index for the total 
system, as well as the partial order indices for some of the components. For example, it is often 
of interest the investigation of the properties of a particular kind of atoms entering a complex 
molecule [30] . Then the order indices can be introduced for the studied atoms characterizing the 
level of their ordering inside the molecule. 

The order indices serve as a measure of internal ordering and correlations and could be used 
as an additional characteristic for finite quantum systems. 
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Figure Captions 



Fig. 1. Dimensionless sound velocity s, quasi- condensate fraction no, and anomalous average 
a as functions of the gas parameter 7. 

Fig. 2. Order index uj{pi) as a function of the gas parameter 7 for different number of atoms: 
N = 10 (sohd hue), N = 10^ (dashed hue), N = 10^ (dotted line), and N = W (dashed-dotted 
line) . 

Fig. 3. Order index u}{pi) as a function of the number of atoms for different gas parameters: 
7 = 0.001 (upper solid line), 7 = 0.1 (dashed line), 7 = 0.3 (dotted line), 7 = 0.5 (dashed-dotted 
line), and 7 = 1 (lower sohd line). 

Fig. 4. Order index uj{pi) for 7 = 1 as a function of the number of atoms A^, demonstrating 
that this index increases as N grows larger than 200. Numerical calculations show that the index 
tends to unity as A?" — >■ 00. 
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Figure 1: Dimensionless sound velocity s, quasi-condensate fraction no, and anomalous average 
a as functions of the gas parameter 7. 




Figure 2: Order index u}{pi) as a function of the gas parameter 7 for different number of atoms: 
N = 10 (sohd line), N = 10^ (dashed line), N = 10^ (dotted line), and = 10^ (dashed-dotted 
line). 
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Figure 3: Order index uj{pi) as a function of the number of atoms for different gas parameters: 
7 = 0.001 (upper solid line), 7 = 0.1 (dashed line), 7 = 0.3 (dotted line), 7 = 0.5 (dashed-dotted 
line), and 7 = 1 (lower solid line). 
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Figure 4: Order index w(pi) for 7 = 1 as a function of the number of atoms A^, demonstrating 
that this index increases as grows larger than 200. Numerical calculations show that the index 
tends to unity as — )■ 00. 
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